Second-harmonic generation in vortex- induced waveguides 
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We study the second-harmonic generation and locahzation of light in a reconfigurable waveguide 
induced by an optical vortex soliton in a defocusing Kerr medium. We show that the vortex-induced 
waveguide greatly improves conversion efficiency from the fundamental to the second harmonic field. 
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Spatial optical solitons have become a topic of active 
research promising many realistic applications and open- 
ing new directions in nonlinear physics I*] . In its simplest 
form, a spatial soliton is a single self-guided beam of a 
specific polarization and frequency. Two (or more) mu- 
tually trapped components with different polarizations 
or frequencies can forma vector soliton. 

One important application of spatial optical solitons 
is to induce stable nondiffractive steerable waveguides 
that can guide and direct another beam, thus creating 
a reconfigurable all-optical circuit. The soliton-induced 
optical waveguides have been studied theoretically and 
demonstrated experimentally in many settings 
It was also shown that soliton waveguides can be used 
for a number of important applications, including the 
second-harmonic generation [g, directional couplers and 
beam splitters 7|, and optical parametric oscillators 0. 

Dark solitons - localized dips on a background inten- 
sity - are more attractive for soliton waveguiding applica- 
tions because of their greater stability and steerability T| . 
Their two-dimensional generalization, optical vortex soli- 
tons 1^, may have a number of potential advantages. Op- 
tical vortex solitons are light beams self-trapped in two 
spatial dimensions and carrying a phase dislocation. A 
systematic analysis of the w aveg uides created by vortex 
solitons in a Kerr medium [ifl [lil demonstrates that 
an optical vortex can guide both weak and strong probe 
beams, and that in the latter case the vortex creates a 
stable vector soliton with its guided component . 

In this Letter, we study the second-harmonic gener- 
ation in a reconfigurable vortex-induced waveguide and 
determine conditions for significant enhancement of the 
conversion efficiency. We also describe novel types of 
three- component vector solitons created by a vortex beam 
together with both fundamental and second-harmonic 
parametrically coupled localized modes guided by the 
vortex-induced waveguide. 

We consider two incoherently coupled beams with fre- 
quencies Wo and uji propagating in a bulk nonlinear Kerr 
medium. The wg-beam propagates in a self-defocusing 
regime and carries a phase dislocation. We assume that 
the phase-matching conditions of the second-harmonic 
generation (SHG) are fulfilled for the fundamental wave 
of frequency uJi guided by the vortex waveguide, so that 
it generates a second-harmonic (SH) wave with the fre- 
quency 2u;i. The SH wave is parametrically coupled 



to the fundamental one and is also guided by the vor- 
tex waveguide. Evolution of the slowly varying beam 
envelopes of the vortex beam, the fundamental guided 
wave, and the SH wave can be described by the following 
system of three coupled dimensionless equations 



i^r + - + c^kP + P\v\'^)u = 0, 
oz 

I— — h /\±w -\- w V — a\u\ w = v, 
oz 

du 1 
2i— + A^i; - Pv + -w^ - = 0. 



(1) 



where it, w, and v are the normalized slowly varying 
complex envelopes of the vortex beam, the fundamen- 
tal field, and the SH field, respectively. Other nota- 
tions are: the Laplacian Aj^ refers to the transverse 
coordinate r = {x,y) measured in units of rn, where 
rl = 3x^^Vl6w2[x(2)]2 (see details in Ref. 12)), z is 
the beam propagation coordinate measured in units of 
zq — 2kirQ. The parameter /S = 2zoAk is proportional 
to the wavevector mismatch Afc = 2ki — k2, whereas the 
nonlinear coupling coefficients a and p are proportional 
to the corresponding third-order tensor components jj], 
and the self-action effects for the fundamental and SH 
fields are neglected. Equations ^ are valid when spatial 
walk-off is negligible and the fundamental frequency coi 
and its second harmonic are far from resonance. 

We emphasize that the model ^ is the simplest of its 
kind, which is most suitable for our feasibility study of 
SHG in vortex-induced waveguides. It is clear that mod- 
elling of particular experimental setups for realization of 
this concept would require modifications of Eqs. (^3), ac- 
cording to the geometry of an experiment and properties 
of nonlinear materials. For example, in photorefractive 
crystals @ one should take into account the nonlinearity 
saturation effect. 

First, we analyze the stationary solutions of the model 
in the form of the (2+l)-dimensional radially sym- 
metric nonlinear modes. We look for spatially local- 
ized solutions in the polar coordinates (r, cf)) of the form 
u = ■u(r)e~*^e*'^, w — w{r)e^^^ , and v — v{r)e^'^^^, with 
the following asymptotic: u{r) 1, and (v{r), w{r)) 
for r 
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oo. Then, the mode amplitudes sat- 



2 




CD 
■D 

"q. 





p = 1 p 


_ w\ 









10 



CD 
T3 
3 



J \w 


p = -1 

;l = 1.0 


E 






u 










FIG. 1: Spatial profiles of the three- wave vector soliton com- 
ponents for the points A to F marked in Fig. H Shown are: 
the vortex amplitude u{r) (thin solid), the fundamental field 
w{r) (thick solid), and the SH field v{r) (dashed) at the indi- 
cated values of /? and A. 



isfy the system of z-independent equations 



1 



ArU ttu + m - (u-^ + 2w^ + 8v^)u = 0, 

ArW — Xw + wv — 2u^w = 0, 

ArV - (4A + (3)v + ^w^ - 8u^v = 0. 



(2) 



where = {l/r)d/dr{rd/dr) is the radial part of the 
Laplacian, and for definiteness we have specified the 
parameters of the cross-phase modulation interaction, 
tJ = 2 and p — 8. In Eqs. Q, the real propagation con- 
stant A must be above cutoff, A > Ac = max(0, — /3/4), 
for w and v to be exponentially localized. 

Using the standard relaxation numerical technique, we 
find the families of localized solutions of the system 
for allowed values of /3 and A. In Fig. ^ we show several 
examples of the profiles of the three- wave localized solu- 
tions for selected values of the parameters /3 and A. The 
numerical results are summarized in Fig. |21 which shows 
the existence domain as a shaded region of the plane 
(A, /3) with the boundary A — Xth (solid curve) found nu- 
merically and the asymptotic lines A = and A = —(3/4: 
(dashed) found from a simple analysis of Eqs. 

All three- wave solutions of Eqs. |(5J can formally be 
divided into two categories according to the dominant 
regime of their formation: (i) vortex- waveguiding regime, 
A < 0, and (ii) quadratic solitons regime, A > . For 
Ath < A < 0, the parametrically coupled modes w and v 
are localized only in the presence of the vortex, and can 
be regarded as two guided modes of the effective vortex- 
induced waveguide. Examples of such solutions are pre- 
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FIG. 2: Region of existence (shaded) of the three-component 
vector solitons of the model Q in the plane (A,/?). Marked 
points correspond to the localized modes shown in Fig. 



sented in Fig.^by the cases A, C, and D. For A > 0, the 
positive values of the propagation constant correspond to 
an effectively self-localizing parametric nonlinearity act- 
ing between the fields w and v. These components can 
then become localized even without the vortex, in the 
form of a parametric quadratic soliton (see Fig. ^ the 
cases B, E, and F). Parametric coupling between the two 
fields is defined, as expected, by the value of the phase- 
matching parameter /3. 

In order to study the SHG process in the vortex- 
induced waveguide, we employ the stationary solutions 
obtained above and analyze numerically the evolution of 
the beams in the case when the SH component is ab- 
sent at the input. We perform all our calculations for 
the case of a finite- extent input vortex beam, obtained 
by superimposing the stationary vortex profile u{r) onto 
a broad super-Gaussian beam, Usg — u{r) exp[—{r^ /d], 
where d — 10*. This form of initial conditions makes our 
predictions more suitable for experimental verifications. 

The numerical results indicate that the generation of 
the SH field from such an input differs dramatically for 
the vortex- waveguiding and quadratic soliton regimes. 
Indeed, for A < we observe a good correspondence with 
the SHG theory. For large /3, the generated SH field is 
weak and the process corresponds to the so-called nonde- 
pleted pump approximation in the SHG theory. Almost 
perfect SHG is observed for /? close to zero, and in all 
such cases the distortion of the vortex waveguide is weak. 
Figure |3f upper row) and Fig. 0] show an example of the 
SHG process with u{r) corresponding to the point D in 
Fig. 12 A good confinement of both fundamental and SH 
guided modes can be seen with a very good conversion 
efficiency and weak distortion of the vortex beam. 

However, in the quadratic soliton regime, when A > 0, 
the strong parametric interaction between the guided 
components does not allow good energy conversion be- 
tween the harmonics. Instead, even for a high-intensity 
fundamental input, both the fundamental and SH fields 
approach a stationary state with nonzero but low- 
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FIG. 3: Examples of SHG in the vortex-induced waveguide 
with no SH field at the input and the parameters correspond- 
ing to the point D (upper row) and point E (lower row) in 
Fig.0and|5| Notice the scale differences between the top and 
bottom rows. 
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FIG. 4: Grey-scaled images of the vortex waveguide and the 
guided modes for the SHG process. Initial conditions corre- 
spond to a vortex carried by a Gaussian beam and the fun- 
damental wave, both corresponding to the point D in Fig. |5| 

amplitude components. The SHG process becomes even 
worse for the negative phase-matching. Figure (^flower 
row) shows an example of a very strong mode coupling 
and vortex distortion corresponding to the parameter re- 



gion /3 < and A > 0. 

If the vortex is removed at the input in the vortex- 
waveguiding regime (A < 0), the SHG conversion effi- 
ciency drops by at least one order of magnitude or more, 
and both the strong fundamental and weak SH fields 
diffract rapidly. In the quadratic soliton regime (A > 0), 
the effective self-focusing nonlinearity of the second-order 
parametric interaction between the fields w and v allows 
the formation of two- wave parametric solitons even with- 
out the vortex component. However, in this case the 
input power does not transfer into the SH field, it un- 
dergoes a redistribution between the harmonics in such 
a way that both fields either approach a stationary state 
corresponding to a (2-|-l)-dimensional quadratic soliton 
(above the existence threshold), or just diffract (below 
the threshold). To summarize, our study of SHG in 
vortex-induced waveguides in different regimes suggests 
that the enhanced conversion efficiency can be achieved 
only in the vortex-waveguiding regime. 

Possible experimental realizations of the concept of the 
SHG in vortex-induced waveguides can be achieved in a 
crystal of FeiLiNbOa where phase-matching can be satis- 
fied through the birefringence effect at the angle = 81° 
with respect to the z-axis, provided the four-wave mix- 
ing effect is suppressed. The other possibility is to em- 
ploy photorefractive crystals and the temperature tuning 
technique, similar to that reported earlier Q. 

In conclusion, we have analyzed the simultaneous 
guidance of both the fundamental and second-harmonic 
waves by an optical vortex soliton. We have described 
novel classes of three-wave parametric solitons with a 
vortex-soliton component, and have studied the second- 
harmonic generation in the vortex-induced waveguides. 
For the first time to our knowledge, we demonstrated 
that larger conversion efficiency of the SHG process can 
be achieved in the vortex-waveguiding regime. 
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